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Crystallization of trions in SU(3) cold atom gases trapped in optical lattices 
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Cold Femiionic atoms with three different hyperfine states confined in optical lattices show pro- 
nounced Atomic Density Waves (ADWs). These ADWs are pinned due to the confining potential 
that traps the atoms in the optical lattice and can be considered a crystal of strongly bound trions. 
We show that the crystalline phase is incompressible and robust against SU(3) symmetry breaking 
interaction. We also show that it is generic in the presence of the trap due to its incompressible 
characteristics. A simulation of the time evolution of the fermionic gas after the trapping potential 
is switch off shows that the formation of the trionic crystal is marked by a freezing of the dynamics 
in the lattice. 
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Ultracold Bose or Fermi gases confined in artificial 
optical lattices have become an important instrument 
for investigating the physics of strong correlations. The 
parameters and dimensionality of these systems can be 
tuned with very high precision and unprecedent control 
[l| by means of an atomic Feshbach resonance or in an 
optical lattice by changing the depth of the wells. One of 
the main achievements in this subject is the experimental 
observation of a superfluid to Mott insulator transition 
in a three-dimensional optical lattice with bosonic *^i?6 
atoms y . Very interesting experimental results have also 
been obtained for fermions [3|]. 

High spin fermions can be directly studied with cold 
atoms in more than two hyperfine states. This kind of 
systems could give rise to new exotic phases in optically 
trapped cold atoms. The general structure of high-spin 
Cooper pairs has already been theoretically addressed 
[J| . More recently, the emergence of triplets and quartets 
(three or four fermion bound states) has been explored 
[5|, [y, III, [SJ, |9|. In the case of SU(3), fermions with three 
hyperfine (color) states, these investigations have sug- 
gested that it would be possible to observe the transition 
from baryonic matter at low densities to a color super- 
conducting state at higher densities [8|, |lO| . At least two 
alkali atoms ^Li and '^'^K seem to be possible candidates 
for the experimental realization of a SU(3) fermionic lat- 
tice with attractive interactions [a]. 

The paradigmatic model used for exploring these 
exotic features is the one-dimensional N-component 
fermionic Hubbard model with a local attractive inter- 
action [lO| . A low-energy effective theory can be derived 
through linearization of the dispersion relation like in the 
usual SU(2) Hubbard Hamiltonian ll|, [ij]- This theory 



predicts new scattering channels in the half- filled case, 
i.e. the Umklapp scattering terms that do not conserve 
momentum but conserve quasi-momcntum. In a recent 
work by Capponi et al. 113] the bosonization parameters 
of SU(N) theories have been investigated by means of 



large scale DMRG calculations. Regarding SU(3) their 
conclusion was that, for low densities, the dominant in- 
stability is that of a trionic superfluid and that color 
superconductivity is strongly suppressed. There is a 
competition between the superfluid trionic behavior and 
ADW order with the former phase dominating for generic 
fillings, and the latter probably appearing at half-filling. 

The main goal of this letter is to study the emergence 
of a trionic crystalline phase in three color atoms loaded 
in optical lattices. We will show how the harmonic con- 
finement potential can pin down an ADW at high den- 
sities in optical lattices, that is the precursor of crystal 
phase of trions (CPT). As the system gains energy due 
to local Umklapp scattering, the confined fermionic cloud 
deforms itself in the lattice to be at half-filling in the cen- 
ter of the trap. Similar effects have also been reported 
in the repulsive and attractive SU(2) cases [lj,ll5|, how- 
ever for three color systems the appearance of an ADW 
is enhanced dramatically. The critical value of the at- 
tractive interaction U for the CPT to be fully formed 
is approximately U — —4. Two main properties of this 
phase appear: a) it is an incompressible state, and b) the 
dynamics is almost frozen due to the suppression of the 
hopping of the trions between lattice wells. Fermionic 
color pairs are only dominant at weak coupling (a small 
region around U = —1). 

The low energy physics of cold fermionic atoms with 
three different hyperfine states trapped in an optical lat- 
tice can be described by an SU(3) generalization of the 
Hubbard Hamiltonian, 



H = -tJ2 (/U" + h-c.) + E 
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The sums over a and j3 go over the three colors The 
operators //^ and fia are the creation and destruction 
operators of an atom in site i with color a. We consider 





FIG. 1: (Color online) Density for different values of U for 
L = 60, iV = 36, and V = 0.001333. 



different values of the on-site interaction between the dif- 
ferent color pairs Uaj3 to be able to include SU(3) sym- 
metry breaking terms as they would probably be present 
in future experiments. The site label i goes from to 
L — 1, with L being the total number of lattice sites. All 
energies are expressed in units of i (t = 1). 

We study the ground state properties of the Hamil- 
tonian ^ with the DMRG algorithm [l6| that provides 
very accurate numerical results. The modifications of the 
DMRG procedure due to the trap potential that breaks 
the translational symmetry have been described in pre- 
vious publications [la . llTI . Il8l | . In order to obtain enough 
accuracy we have kept up to 1200 states in each iteration. 

In Figure [T] we show the atom density for equal color 
populations N = Nr + Ng + Nb = 36, parameters 
V = 0.003 , L = 60, and different values of the attrac- 
tive SU(3) interaction U = Urg = Urb = Ugb- The figure 
clearly shows the formation of the CPT as \U\ increases. 
For [/ < — 4 the phase is fully developed. Increasing \U\ 
compresses the density as interaction energy gains over 
kinetic energy and the particles prefer to be in the center 
of the trap. However, the crystal phase is incompressible 
and, once it is formed, increasing the value of the attrac- 
tion does not compress the particle density any further. 

The formation of these phases can be studied by means 
of local correlation functions that measure the number of 
pairs and the number of trions. For the number of pairs 
we define 



PiaP = {niariif}) - {ntrnigHib) 



(2) 



In the non-polarized we will be interested in the total 
number of pairs at each site i. 
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The local correlation that counts the number of trions is 

Tj = {nirTiigna) . (4) 



FIG. 2: (color online) Extrapolated total number of fermions 
in pairs (2P) and in trions(3T) compared with the single 
fermions(F) for different sizes as a function of \U\. Inset: 
Fitting second order polynomial for 2P/N as a function of 
1/iV for U = -2. 



For non-interacting SU(3) fermions Piaf} — {nia) {nip) — 
{riir) {riig) (nib) and T^ = (n^r) {uig) (nib). 

These local correlations can be summed over all the 
sites on the system to obtain the total number of fermions 
in pairs (2P) and the total number of fermions in trions 
(3T). The number of uncorrelated fermions (F) is then 
F = Q - 2P - 3T. 

The proper finite size scaling of the systems, as de- 
scribed in '215, is given by N'^V = constant when 
N -^ 00 and V^ ^ 0. We checked the validity this scaling 
law by comparing results for different sizes. For a given 
value of U the results for the fraction of fermions in tri- 
ons 3T/N and the fraction of fermions in pairs 2P/N 
as a function of the inverse size 1/A^ can be easily fit- 
ted by a second order polynomial. In this way we can 
extrapolate the values of trion and pair fractions in the 
thermodynamic limit. Figure [2] shows these results. It 
can be seen that the number of trions increases rapidly 
with the interaction strength until almost all fermions are 
present as trions for [/ < — 4. The number of pairs is only 
dominant in a small region around U = —1 suggesting a 
color superconducting phase. An example of the quality 
of the extrapolation procedure is shown in the inset of 
the figure. 

To understand better the origin of the ADW state we 
will compare the DMRG ground state energy with sev- 
eral approximations. In the strong interaction limit it 
is possible to treat the hopping term as a perturbation. 
Ignoring the effect of the trap, the zeroth order energy of 
each trion is 3f7. The first non-zero correction appears 
in second order perturbation theory giving a lowering 
of 3t^ /U for two next to nearest-neighbor trions, due to 
quantum fluctuations of the fermions within the trions. If 
two trions sit in neighboring sites this term is completely 
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FIG. 3: (Color online) Absolute value of the ground state 
energy |_Eo| in logarithmic scale as a function of |C/| for the 
case with L = 40, 1/ = 0.003, and iV = 24. DMRG results 
are compared with a HF approximation, a CPT, and a Fermi 
liquid state of trapped trions. 



suppressed. Therefore, it acts as an effective repulsion 
between nearest-neighbor trions, being at the origin of 
the CPT at half-fiUing [19^. UnHke the SU(2) case, the 
first contribution to hopping appears in the third order 
of perturbations, thus a CPT will always dominate for 
strong enough attraction. 

Adding the effect of the trap we have calculated the 



energy of the CPT state, defined as Hi' ' 






0), 



where the i' labels the even sites closer to the trap cen- 
ter until all the fermions are accounted for. We have 
also calculated the energy of a Fermi liquid (FL) state, 
defined as a Slater determinant in the basis that diag- 
onalizes the hopping term plus the trapping potential. 
The CPT state valid in the strong coupling limit and the 
FL state appropriate for weak coupling limit are both 
Slater determinants independent of the attractive inter- 
action U . Therefore, it is expected that a Hartree-Fock 
(HF) approximation that looks for the best Slater deter- 
minant for a given value of U will interpolate between 
the two limits, and it will provide a better approxima- 
tion for intermediate interactions. In Fig. [3] the ground 
state energy i?o as a function of \U\ for the three differ- 
ent approximations is compared with the DMRG results. 
As expected, the HF energies are very accurate for small 
and large values of \U\. They deviate from the DMRG 
results appreciably in a region around f7 = — 2 where 
the intermediate color superconductivity phase appears 
to take place. 

When the effect of the trap is not taken into ac- 
count the ADW order is non-generic occurring only at 
exactly half-filling. According to our simulations the 
CPT state is generic in the presence of a trap and oc- 
curs for a wide range of values of the number of particles 
in the trap, the confining potential, and the interaction 



FIG. 4: (Color online) Normalized size of the fermionic cloud 
an/N as a function of N'^V for different values of U and A^ = 
36. 



strength. To address this issue we have calculated the 
extension of the fermionic cloud in the trap, defined as 
cr„ = \/Si "•«(* ^ ^/2)^: as we vary the confining po- 
tential for different values of the interaction strength. 
Results of these calculations for A^ = 36 are shown in 
Fig. m The horizontal dotted line denotes the value of 
cr„ for the CPT while the continuous horizontal line at 
the bottom denotes the Band Insulator (BI) state with 
all atoms filling the bottom of the trap. For U < —4, 
when the value of the CPT is reached, there is a plateau 
of CT„ as a function of N'^V marking the appearance of 
an incompressible state. For U < —6, new plateaus for 
lower (T„ and higher values of V appear. These other 
plateaus correspond to the presence of a band insulator 
in the center and symmetric CPTs at the sides of the 
trap. The crossover from the CPT to the BI as we in- 
crease the confinement consists on a number of steps in 
which each of the trions in the CPTs in the outer side 
of the cloud is forced to move the BI in the center (an 
analogous phenomenom has been recently described for 
spin 3/2 atoms |9|]). 

In what follows we will explore the evolution of the 
CPT when the trap is opened. We expect that its dynam- 
ics should be frozen since the effective hopping of a bound 
trion is extremely small (third order in perturbation the- 
ory). We have prepared a state of five trions made out of 
15 fermions in a lattice of size L = 40 and v = 0.005. The 
harmonic confinement potential is switched off at t = 0, 
and the state evolves with the rest of the Hamiltonian 
(HI). The temporal evolution is computed with the time 
dependent DMRG algorithm [20| . In Fig. [5] we show the 
evolution in natural time units of the density for two se- 
lected values of [/: U — —2 in the color superconducting 
phase and U = —8 well in the CPT. As a general trend, 
we observe that the expansion of the cloud almost ceases 
for U < —6. In order to quantify this statement we have 




FIG. 5: (Color online) Evolution of the ground state density 
after a sudden switch off of the confinement potential for the 
case with L = 40, V = 0.005, and Q — 15. Left panel: 
U = -2. Right panel: U = -8. 
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FIG. 6: (Color online) Velocity of expansion as a function of 
the interaction in a case with A'^ = 15. 



calculated the velocity of the most external trion as a 
function of U . The velocity is computed by means of a 
hnear fit of the position of the external peak as a function 
of time. The results are depicted in Fig. [6] where can see 
how the velocity is strongly suppressed for U < —4. 

As already mentioned, a possible candidate for real- 
izing the trionic phase is ^Li for which the magnetic 
field dependence of the three scattering lengths has been 
measured [22l |. The attractive U interaction in some ap- 
propriate regions of magnetic field can be estimated as 
Urg = Uq, Urb = 1.23C/0, and Ugb = IMUq d. Further- 
more, the effective value of Uq can be controlled changing 
the depth of the optical lattice wells. We have numeri- 
cally checked the effects of this experimental anisotropy 
in the interaction. The CPT persists without noticeable 
changes even for stronger SU(3) symmetry breaking in- 
teractions. 

In summary, we have shown that an CPT of trions 
is a very robust and generic phase of atoms with three 
hyperfine states confined in optical lattices. Due to the 
incompressible properties of this state it should be easily 
observed for experiments with high-density of fermionic 



atoms. The appearance of this state is marked by a freez- 
ing of the size of the cloud when the trapping potential 
is switched off, as we have shown through a full time 
simulations of the evolution. In future experiments of ul- 
tracold fermionic atoms the existence of the CPT could 
be revealed by light-scattering diffraction experiments. 
Moreover, radio- frequency spectroscopy could be used to 
probe the dissociation dynamics of the trions. 
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